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Abstract 

Elastic and spin-changing inelastic collision cross sections are presented for cold and ultracold 
magnetically trapped NH. The cross sections are obtained from coupled-channel scattering cal- 
culations as a function of energy and magnetic field. We specifically investigate the influence of 
the intramolecular spin-spin, spin-rotation, and intermolecular magnetic dipole coupling on the 
collision dynamics. It is shown that ^^NH is a very suitable candidate for evaporative cooling 
experiments. The dominant trap-loss mechanism in the ultracold regime originates from the inter- 
molecular dipolar coupling term. At higher energies and fields, intramolecular spin-spin coupling 
becomes increasingly important. Our qualitative results and conclusions are fairly independent of 
the exact form of the potential and of the size of the channel basis set. 
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I. INTRODUCTION 



Cold (T < 1 K) and ultracold (T < 1 mK) molecules have attracted great interest in 
recent years due to their potential apphcations in condensed matter physics II, quantum 
computing j^, 3], high-precision spectroscopy and physical chemistry [7-ll|. Experi- 



mental methods for producing ultracold molecular samples include photoassociation 12] and 
Feshbach association in which molecules are formed by linking two ultracold atoms. 
These indirect cooling techniques are, however, limited to molecules consisting of two (alkali 
metal) atoms that can be efficientl y p re-cooled to ultralow temperatures. Direct cooling 
methods such as buffer-gas cooling IJ] and Stark deceleration 15|], in which molecules are 
cooled from room temperature, are applicable to a much wider range of species. 

The NH radical is an excellent candidate for direct cooling experiments, due to its rel- 
atively large magnetic and electric dipole moments. Ground-state NH(X^S~) has been 



successfully cooled with a helium buffer gas and confined in a magnetic trap 



Id, 



16- 



8l. Metastable NH(a^A) has been Stark decelerated and trapped in an electrostatic field 
Furthermore, an optical pumping scheme has been developed to transfer the Stark- 
decelerated molecules to the electronic ground state and subsequently accumulate them in 
a superimposed magnetic trap [20|. 

At present, direct cooling experiments for NH have achieved temperatures of a few hun- 
dred mK. One of the key challenges is to produce molecules in the ultracold regime, which 
requires a second-stage cooling method. Examples of such methods include sympathetic 



cooling with an ultracold co-trapped species and evaporative cooling. The latter 

approach relies on elastic NH -|- NH collisions that thermalize the gas as the magnetic trap 
depth is slowly decreased. Inelastic collisions, however, can induce Zeeman relaxation and 
consequently lead to trap loss. Therefore, in order for evaporative cooling to work, elastic 
collisions should be much more efficient than inelastic ones, typically by a few orders of mag- 
nitude 2^1. Evaporative cooling has already been achieved for buffer-gas cooled metastable 



helium 



271], indicating that the densities in buffer-gas cooling experiments are high enough 



to obtain a sufficiently high elastic collision rate. This also provides hope for other buffer-gas 
cooled species such as NH. 

A collision complex consisting of two magnetically trapped NH molecules is in the 
\S = 2, Ms = 2) state, with S denoting the total spin and Ms its laboratory-frame pro- 
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jection. Inelastic NH + NH collisions can change the Ms quantum number and/or the total 
spin of the dimer. If the total spin is flipped, the complex undergoes a transition from the 



quintet state to the chemically reactive singlet or triplet states |28l-l30| . Although unfavor- 



able for evaporative coo 



controlled chemistry 



ing, the latter two spin states may be relevant in the context of cold 



Previous theoretical work by Kajita 3l|] and Janssen et al. 32| has shown that NH 
is a promising candidate for molecular evaporative cooling, particularly the bosonic ^^NH 
isotope. Fermionic ^^NH may also be cooled into the /xK regime either by applying an 



electric field 



3l| or by preparing the molecules in different hyperfine states 32|. It must 



be noted, however, that the results of Kajita were obtained from approximate analytical 
methods in which only the electric dipole-dipole and magnetic dipole-dipole interactions 
were included. The work of Ref. 32|] involved rigorous quantum scattering calculations on a 



full 4-dimensional potential-energy surface, but focused only on NH + NH collisions in zero 
field. 

In this paper, we consider collisions between ^^NH molecules in the presence of an ex- 
ternal magnetic field. We present full quantum scattering calculations on an accurate ah 
initio quintet potential to investigate the collision dynamics at low and ultralow tempera- 
tures. Intramolecular spin-spin coupling, spin-rotation, and intermolecular magnetic dipolar 
coupling are explicitly included in the calculations. We will identify the dominant trap-loss 
mechanisms and provide a detailed discussion of the dynamics. 



II. CALCULATIONS 



A. Theory 

We consider elastic and inelastic collisions between two identical NH(X^S~) molecules 
in the presence of a magnetic field. Since bosonic ^^NH is more suitable for evaporative 
cooling that fermionic ^^NH, we will focus only on the ^^NH isotopologue. We treat the 
colliding molecules as rigid rotors. The coordinate system is defined in a space-fixed frame, 
with R = {R, Q) and Q = ($, 0) denoting the intermolecular vector between the centers 
of mass of the molecules. The coordinates loa = {Oa,(Pa) and ub = {(^b,<Pb) are the polar 
angles of monomers A and B, respectively. We neglect hyperfine structure and assume that 
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both molecules are in their nuclear-spin stretched states |/, M/ = /), with J = /n + /r = 1 
denoting the maximum total nuclear spin and Mi its laboratory-frame projection. The 
scattering Hamiltonian is then given by 

fi^ d"^ L"^ .... 

^ ^ ~2jiR'dR^^ ~^ 2/IR2 ~^ ^s{R,(^A,(^b) + Knagn.dip(^, 5'^, Sb) + Ha + Hb, (1) 

where is the reduced mass of the complex, is the angular momentum operator associated 
with end-over-end rotation of the vector R, Vs{R,U!a,'^b) is the potential-energy surface 
for total spin S, Knagn.dip(-R, Sa, Sb) is the intermolecular magnetic dipole interaction, and 
Ha and Hb are the Hamiltonians of the two monomers. The intermolecular magnetic dipole 
term is given by 

2 

Wdip(J2, Sa, Sb) = -V^glf^l^ 5^(-l)«C2,-,(f^)[^^ ® ^^if (2) 

where gs ~ 2.0023 is the electron g-factor, /i^ is the Bohr magneton, a is the fine-structure 
constant, and the factor in square brackets is the tensor product of the monomer spin 
operators Sa and Sb- The monomer operators Hi {i = A,B) each contain a rotation, 
spin-rotation, intramolecular spin-spin, and Zeeman term, 

H, = BoN^ + 7iV. ■ S, + ^VeAss 5^(-l)«C2,-,(u;.)['5^ ® ) + g^^^B ■ S„ (3) 

where Ni is the operator for the rotational angular momentum of monomer i and B is the 
magnetic field vector. The numerical values for the rotational, spin-rotation, and spin-spin 
constants of ^^NH can be obtained from those of ^"^NH 33| using isotope scaling (see e.g. p. 



239 of Ref. |3|): Bq = 16.270 340 cm'\ 7 = -0.054 60 cm'^, and Ags = 0.919 89 cm-\ 

We will assume that all three spin states of the complex are described by the nonreactive 
quintet surface, i.e. Vs(R, uja,i-or) = V2{R,coa,ojb)- Field-free calculations have shown that 



this assumption is reasonable 



32| . The 5 = 2 surface is taken from Ref. (29|. The potential 



is expanded in coupled spherical harmonics Yl,m j35| . 

V{R,Ua,U;b)= Y1 VLA,LB,LABiR)^LA,LB,LABi^,^A,UJB), (4) 
La,Lb,Lab 

ALA,LB,LABi^^^A,UJB) = Yl {LaMaLbMb\LabMab) 

MaMbMab 

X YLA,MAi^A)YLB,MBi^B)Yl^^^MAB(^)^ (5) 
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where {LaMaLbMb\LabMab) is a Clebsch-Gordan coefficient and the superscript * denotes 
complex conjugation. The subscript 5* = 2 has been omitted for brevity. As noted previously 
in Ref. |32|, the t^L^,La,LAs(-R) expansion coefficients of Ref. [29| must be multiplied by a 
factor of (-l)^-4--^s(47r)-3/2(2L^B + 1)[(2L^ + 1)(2Lb + l)]^/'^ to obtain the potential in the 
form of Eq. (jl]) . 



We expand the wave function in a space-fixed uncoupled basis introduced in Ref. 
The angular functions are written as products of the eigenfunctions of Nf, Ni^, Sf, Si^, L^, 
and 1/2, 

\NAMN^)\SAMsJ\NBMr,,)\SBMs,)\LML) = ItaTb)!^^^), (6) 

with Na and Nb ranging from to A^max and L = 0, . . . , Lmax- The matrix elements of the 
scattering Hamiltonian [Eq. ([1])] in the basis of Eq. (|6]) have all been given by Krems and 
Dalgarno 36|. Note that the factor of g^f^'s ~ 1-0023 e^fi^/ml has been neglected in their 



expression for the magnetic dipole interaction. Furthermore, it should be taken into account 
that our potential expansion [Eg (|1D] differs by a factor of {-l)^^'^B{An)^/'^{2LAB + 1)^''^ 



from the expansion used in Ref. 



Since the monomers are identical, we can exploit the permutation symmetry of the sys- 



tem to minimize the computational cost. Following Ref. 37|], we employ a normalized, 
symmetrized basis of the form 

\K,,lmJ = ^ [hAlB)\LM,) + r^(-l)^|7B7^)I^M^)] , (7) 

with ?7 = ±1 defining the symmetry of the wave function with respect to molecular inter- 
change. It is sufficient to restrict the basis to a well-ordered set of molecular states, i.e. 
7a > 7_B- Finally, we note that the parity, e = (—1)^-4+^^3+-^^ and the total angular momen- 
tum projection quantum number, Ai = Mjy^ + Ms^ + M^g + Msg + M^, are also conserved 
during collision. Thus, the scattering calculations may be performed for a single value of 
Ai, rj, and e. Here we consider only collisions between bosonic ^^NH molecules in their 
spin-stretched and rotational ground state (A^^ = Nb = 0), for which we have rj = +1 and 
e = +1. Note that the ffist excited rotational state has an energy of ~ 32 cm~^ (46 K), and 
is therefore inaccessible at the energies considered in this work. 



We solve the cou 
and Manolopoulos [38 



ed equations using the hybrid log-derivative method of Alexander 
, which uses a fixed-step-size log-derivative propagator in the short 



range and a variable-step-size Airy propagator at long range. Matching to asymptotic 



boundary conditions yields the scattering 5'-matrices, from which the cross sections can 
be readily obtained. We note that, due to the intramolecular spin- rotation and spin-spin 
couplings, the basis functions of Eq. ([7]) are not exact eigenf unctions of the asymptotic 
Hamiltonian Ha + -^b, while the 5*- matrices must be constructed in terms of these eigen- 
functions. As detailed in Refs. 36| and 37|, an additional basis transformation of the 



log-derivative matrix is therefore required before matching to the asymptotic boundary 
conditions. We will denote the symmetrized channel eigenf unctions as with 
|7^7b) = \{Na, Sa)Ja, Mj^)\{Nb, Sb)Jb, Mjg) defining the molecular eigenstates. Here iVj 
is used instead of A^j, because Ni is strictly not an exact quantum number due to the in- 
tramolecular spin-spin coupling. However, the intramolecular coupling is relatively weak and 
Ni,Mi^., and Ms- may be treated as approximately good quantum numbers. Specifically, 
for the ^^NH rotational ground state, the magnetically trapped | Jj = 1, Mj. = 1) component 
contains 99.992% of | Aj = 0, Mjy- = 0, Si = 1, M^- = 1) for all fields considered in this work. 
The cross sections at total energy E are calculated using the expression ^\ 

•^lA^B l^Ml L',M' 



where ky^y^ is the wavenumber for the incident channel, k'^^y^ = 2fi(E — e^^ — eyg)/h'^, 
e^^ + e^g is the corresponding channel energy, and the T-matrix elements are defined as 

B. Computational details 

We performed the scattering calculations for different magnetic fields using the 



MOLSCAT package 



39 



40|, extended to handle the basis set of Eq. ([7]). The log- derivative 
propagation was carried out on a radial grid ranging from 4.5 to 15 in steps of 0.02 
oq. The Airy propagation ranged from 15 to 50 000 oq. We included basis functions up to 
Amax = 2 and Lmax = 6, yielding a maximum number of 1038 channels for a single calcu- 
lation {M. = 0). The expansion of the quintet potential-energy surface was truncated at 
La = Lb = Q. As mentioned in Section [II At the chemically reactive singlet and triplet 
potentials were not included in the calculations. 
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III. RESULTS AND DISCUSSION 



A. Adiabatic potential-energy curves 

Before discussing the cross sections, we first consider the adiabatic potential curves. 
These are obtained by diagonahzing the interaction matrix at fixed R over a grid of R- 
values and subsequently connecting the corresponding eigenvalues. Asymptotically, these 
adiabatic curves correlate to the molecular eigenstates 7^ and 75. Figure [1] shows the 
long range of the lowest adiabats for Ai = 2 and ivmax = 4 at a magnetic field of 0.1 G. We 
present only the curves with exchange symmetry rj = +1 and parity e = +1, for which s-wave 
scattering in the incident spin-stretched channel is allowed. The asymptotic energy splittings 
between the different molecular states originate from the Zeeman term gs^BB{Msj^ +Msg)- 
For M. = 2 and Lmax = 4, the magnetically trapped state with |Mj^ = Mj^ = 1) has 
partial-wave contributions of L = 0, 2, and 4, as can be seen in Fig. [1] The adiabats 
correlating with \Mj^ = l,Mj^ = 0), \Mj^ = l,Mj^ = -1), and \Mj^ = 0,Mj^ = 0) 
have L = 2 and 4 centrifugal barriers, and the adiabats for \Mj^ = 0,Mj^ = —1) and 
\Mj^ = —1, Mj^ = —1) contain only the L = 4 partial wave. For the \Mj^ = 1, Mj^ = — 1) 
and \Mj^ = 0, Mj^ = 0) states, which have identical Zeeman shifts, the degeneracy is 
further lifted by the intramolecular spin-spin coupling. 

It can be seen that several curve crossings occur in the region between R ^ 1500 and 4500 
oq. If we neglect the small energy shifts due to the intramolecular coupling, the crossing 
points Rc are simply given by 

gsf^sBAMs , (9) 

where Li and Lf denote the values of L for the adiabats correlating to the incoming and 
outgoing channels, respectively, and AM5 = M^^ + M^^ — M^P — Mif^ . The corresponding 
energies at which the crossings occur are given by = h^Li{Li + l)/{2fiRl), defined relative 
to the threshold of the incident channel. We find that several of these crossings are narrowly 
avoided due to the presence of the intermolecular magnetic dipole interaction. Inspection 
of Eq. (30) of Ref. or Eq. (A2) of Ref. 411 shows that e.g. the s-wave incident channel 



with \Ms^ = l,Msg = 1,L = 0) is directly coupled with \Ms^ = l,Msg =0,L = 2) and 
\Ms^ = 0,Msg = 0, L = 2) via the magnetic dipole term, giving rise to the corresponding 
avoided crossings. We note, however, that Knagn.dip contains a second-rank tensor in fl and 



can therefore directly couple channels only ii \Li — Lf \ < 2. The spin operators Sa and Sb 
contained in Vmagn.dip are first-rank tensors and, consequently, and Msg may each differ 
at most by 1. Thus, not all crossings are avoided. To take proper account of these avoided 
crossings and account correctly for the magnetic dipole coupling in the ultracold regime, it 
is essential that the radial grid used in the calculations extends beyond the outermost Rc 
value. The influence of the avoided crossings on the collision cross sections will be discussed 
in detail in the following section. 



B. Cross sections 



Figure [2] shows the total elastic and Mj-changing inelastic cross sections for two magnet- 
ically trapped ^^NH molecules (Mj^ = Mjg = 1) as a function of collision energy. The cross 
sections as a function of magnetic field are given in Fig. |3l We find that the elastic cross 
sections are significantly larger than the inelastic ones over a wide range of energies and 
magnetic fields, suggesting that evaporative cooling of ^^NH is likely to be successful. It can 
also be seen that, in the ultracold regime, the total inelastic cross section decreases dramat- 
ically if the magnetic field strength is reduced. Thus, once the cooling process has started in 
the millikelvin regime at relatively high magnetic field, and continues towards lower energies 
as the magnetic trap depth is decreased, the ratio between elastic and Mj-changing cross 
sections will remain very favorable for evaporative cooling to take place. 

In order to identify the main trap-loss mechanism, we have also performed scattering 
calculations where two of the three spin-dependent coupling terms are set to zero. The 
results for Ai = 2 are shown in Figs. H] and [5l We have verified that the Ai = 2 cross 
sections, which allow for s-wave collisions in the incident channel, are dominant at the 
energies and fields considered in this work. It can be seen that the intermolecular spin- 
spin interaction (Knagn.dip) gives the largest contribution to the inelastic cross section over 
a broad range of energies and field strengths, most notably in the ultracold regime. This is 



24 



251]. At higher collision energies 



also the case for cold and ultracold N -|- NH collisions 
and fields, however, the intramolecular coupling terms become increasingly important. In 
particular the intramolecular spin-spin term causes significant trap loss above E ^ 10~^ K 
and B ^ 100 G. Spin-rotation coupling, which vanishes for pure A^^ = Nb = states, has 
only a very small effect on the total cross section. This is consistent with previous work on 
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the He-NH(3S-) system [l6|, ISQ, |42|. 



The importance of the intermolecular magnetic dipole interaction is most easily under- 
stood by considering the adiabatic potential curves. As discussed in Section IIII A[ the 
avoided crossings between the s-wave incoming channel \Mj^ = 1, Mj^ = 1, L = 0) and the 
\Msj^ = 1, Msj^ = 0,L = 2) and \Ms^ = 0, Msg = 0,L = 2) outgoing channels all occur at 
very long range for small to moderate field strengths. Consequently, the spin-flip induced 
by Vmagn.dip cau take place without having to overcome the rf-wave barrier in the outgo- 
ing channels, and hence this process dominates the inelastic cross section in the ultracold 
regime. We also emphasize that, if the avoided-crossing points Rc fall outside the scattering 
propagation region, i.e. if the radial grid is chosen too small, the inelastic cross sections are 
similar to the case where the intermolecular magnetic dipole term is switched off completely. 
This confirms that the spin-flip due to Knagn.dip indeed takes place at long range, or more 
specifically, at R ^ Rc. 

Kajita 3l|| has shown that, in the Born approximation, the inelastic cross section cTj^/ for 
^^NH [Li = 0,Lf = 2) caused by the magnetic dipole interaction should be proportional to 
5^/2 and E~^/^ if the kinetic energy release is relatively large (fcj ^ kf). This is consistent 
with our results obtained from full quantum scattering calculations in the ultracold regime. 
As the collision energy increases, however, the assumption of ki kf breaks down and the 
cross sections deviate from the B^^'^ behaviour. We find numerically that this is the case for 
energies above ~ 10~^ K at nearly all the field strengths considered in this work (see Fig. 
[3]). In a separate publication, we will give a general analytical expression for the inelastic 
cross section due to Vmagn.dip based on the (distorted-wave) Born approximation, and show 
that the numerical and analytical results are in excellent agreement over a wide range of 
energies and fields. 

At collision energies above about 10 mK or high magnetic fields, there is sufficient energy 
to overcome the d-wave centrifugal barrier in the outgoing channels and, as a consequence, 
short-range effects become important. In particular the intramolecular spin-spin coupling 
;erm, which requires strong anisotropy of the potential in order to induce Zeeman relaxation 



16|, 



36j, enhances the inelastic cross section significantly at energies above ~ 10 mK. For 



the intramolecular contributions, we find that the inelastic cross section behaves as S^/^ at 
moderate field strengths and flattens off to a constant value at very small B. Its kinetic 
energy dependence is proportional to E^^/'^ in the ultracold regime and, for small magnetic 
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fields, also shows a region of E'^ behaviour. This result is consistent with the work of 
Volpi and Bohn, who applied the distorted-wave Born approximation to calculate inelastic 



spin-changing collisions induced at short range [43|: 



a.,^j oc ^ gsftsBAMsf^^'/'' . (10) 

Thus, for an s-wave incoming channel (Lj = 0) and d-wave outgoing channel {Lf = 2), 
there is a region of E"^ dependence when the Zeeman shift for the outgoing channel is small 
compared to the collision energy. At very small fields, the magnetic-field dependence is 
negligible and the cross section becomes constant as a function of B. 

We also point out that, in contrast to the intermolecular spin-spin coupling term, 
the intramolecular spin-spin interaction contains second-rank tensors in Sa and Sb 
and therefore directly couples channels where M^^ and each differ by 0, 1, or 

2. Thus, transitions from Ms^ = 1 to M5- = —1 also become allowed in first or- 
der. This is illustrated in Fig. Ej where the state-to-state inelastic cross sections for 
Ai = 2 are plotted as a function of energy. In the ultracold region, which is dom- 
inated by Knagn.dip, ouly the \Ms^ = 1,Msb = 1) -> \Ms^ = l,Msg = 0) and 
\Ms^ = l,Msj^ = 1) — > \Ms^ = 0,Msg = 0) transitions contribute to the inelastic 
cross section. At higher energies, where the intramolecular spin-spin term plays a signifi- 
cant role, transitions from \Ms^ = 1, Msg = I, L = 0) to the \Ms^ = 1, Ms^ = —1, L = 2), 
\Ms^ = 0, Ms^ = — 1, L = 4), and \Ms^ = — 1, Ms^ = — 1, L = 4) channels become increas- 
ingly important. Note that the latter two have g-wave barriers in the exit channel, and 
hence they are strongly suppressed in the low-energy regime. 

In summary, we have established that the dominant trap-loss mechanism for NH in the 
ultracold regime is the intermolecular spin-spin coupling term, which induces Zeeman re- 
laxation at long range. When the kinetic energy in the outgoing channel becomes large, the 
spin-change is also caused by the interplay of the potential anisotropy and the intramolecular 
spin-spin interaction, which acts at short range. 



C. Dependence on potential and basis-set size 

As a final part of our discussion, we consider the sensitivity to the potential and the 
dependence on the size of the channel basis set. It is well known that low-energy collisions 
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are highly sensitive to the exact form of the interaction potential, so that we must carefully 
take into account the uncertainty in the potential-energy surface. We estimate that the 
quintet potential used in the scattering calculations is accurate to within a few percent. 



Following Refs. [24!] and 



32|, we have studied the potential dependence indirectly by eval- 



uating the sensitivity to the reduced mass. The reduced mass was scaled by a factor of A 
(0.9 < A < 1.1), which is essentially equivalent to scaling of the entire potential by A 24 1. 
Figure [7] shows the Ai = 2 cross sections for different scaling factors at a collision energy 
of 10~^ K and a field strength of 1 G. We find that the cross sections exhibit resonance 
structures at certain A values, indicating a high sensitivity to the potential. However, the 
ratio between elastic and inelastic cross sections is much less sensitive to the reduced mass, 
and the prospects for evaporative cooling remain very positive over almost the entire range 
of A. The contributions from the different spin-changing mechanisms, as described in the 
previous section, are also very similar for different reduced masses. Our qualitative results 
and conclusions are thus reasonably independent of the precise form of the potential. 

We use the same scaling approach to investigate the dependence on basis-set size. Figure 
[8] shows the Ai = 2 cross sections as a function of A for different values of A^max and Lmax 
at E = 10"^ K and B = 1 G. The total number of channels in these calculations ranged 
from 901 for A^max = 2, Lmax = 6 up to 2621 for A^max = 3, L^^x = 6. It can be seen that the 



positions of the resonances shift when the basis size is increased, but the 



genera. 



pattern 



32|. We 



remains essentially the same. This is consistent with other work on NH [23|, 
thus conclude that, within the uncertainty limits of A, our quahtative results are not very 
sensitive to the size of the angular basis set. 

Figure |8] also demonstrates that the cross sections are not yet converged with respect to 
A^max and Lmax- In fact, field-free NH-NH calculations suggest that the basis set should 
extend to at least A^max = 6 and L^a.^ = 7 to achieve full convergence |32|, which would 
amount to a maximum of 25 598 channels for = in the present decoupled basis set. Such 
calculations are highly infeasible with the currently available computer power. However, 
taking into account the uncertainty in the potential and the pronounced resonance structure, 
even a fully converged basis set would not give really reliable numerical values. Nonetheless, 
it must be emphasized that unconverged basis sets can give qualitatively rehable results. As 



already discussed in Refs. [32!] and 



231], full basis-set convergence is therefore not strictly 



necessary within the uncertainty limits of the potential. 
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IV. CONCLUSIONS 



We have carried out full quantum scattering calculations to study cold and ultracold 
^^NH - ^^NH collisions in magnetic fields. The elastic and spin-changing cross sections for 
magnetically trapped NH are found to be very favorable for efficient evaporative cooling. We 
have identified the intermolecular spin-spin coupling term as the main trap-loss mechanism 
at low energies and small magnetic fields, while the intramolecular spin-spin term becomes 
increasingly important at higher energies and fields. The effect of spin-rotation coupling is 
almost neghgible. 

Finally, we have demonstrated that the numerical values of the cross sections are very 
sensitive to the details of the potential, but the qualitative results and conclusions are almost 
independent of the exact form of the surface. The size of the angular basis set, which is 
almost impossible to converge for systems such as NH-NH, does not significantly alter the 
results within the uncertainty limits of the potential. This inherent uncertainty in the 
calculated cross sections, however, clearly highlights the need for reliable experimental data. 
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FIG. 1: Lowest adiabatic potential curves for ^^NH - ^^NH, calculated for Lmax = 4 and Ai = 2 
at a magnetic field of 0.1 G. The molecular eigenstates are labeled by |Mj^, Mjg) and refer to the 
well-ordered states {{Na = 0)Ja = 1,Mj^)\{Nb = 0)Jb = l,Mjg). 
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FIG. 2: Elastic and inelastic Mj-changing cross sections for magnetically trapped ^NH as a 
function of collision energy for various magnetic fields. The elastic cross sections are the same for 
all three magnetic field strengths. 
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FIG. 3: Elastic and inelastic Mj-changing cross sections for magnetically trapped ^^NH as a 
function of magnetic field for various collision energies. 
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FIG. 4: Inelastic Mj-changing cross sections for = 2 as a function of collision energy at 
1 G and 100 G. The different curves correspond to calculations with either the intermolecular 
magnetic dipolar coupling ("Knagn.dip")' the intramolecular spin-spin coupling ("Vintra-ss")) or 
the intramolecular spin-rotation coupling ("Mntra-Ns") included in H. 
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FIG. 5: Inelastic Mj-changing cross sections for = 2 as a function of magnetic field at 10 
K and 10"^ K. The different curves are obtained from scattering calculations with either the 
intermolecular magnetic dipolar coupling ("Knagn.dip")i the intramolecular spin-spin coupling 
( Mntra— ss")) or the intramolecular spin-rotation couphng ("Mntra-Ns") switched on. 
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FIG. 7: Elastic and inelastic Mj-changing cross sections (A^ = 2) for magnetically trapped ^^NH 
as a function of the scaling factor A, calculated at a collision energy of 10~^ K and a magnetic field 
strength of 1 G. 
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FIG. 8: Elastic and inelastic Mj-changing cross sections {M = 2) as a function of the scaling 
factor A, calculated for different basis sets at a collision energy of 10^^ K and a magnetic field 
strength of 1 G. 
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